Spin Hamiltonian for which the Chiral Spin Liquid is the Exact Ground State 
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We construct a Hamiltonian that singles out the chiral spin liquid on a square lattice with periodic 
boundary conditions as the exact and, apart from the two-fold topological degeneracy, unique ground 
state. 



Introduction. The fractionalization of quantum num- 
bers, in which the excitations of a strongly-correlated 
system carry only a fraction of the quantum numbers 
of the constituents, is currently of great interest in con- 
densed matter physics and a significant body of recent 
work has focused on finding solvable theoretical models 
in which the phenomenon occurs [J El E|, 0, E| . I n ad- 
dition to its intrinsic interest, the phenomenon of frac- 
tionalization may well have a bearing on one of the most 
vexing problems in condensed matter theory, should the 
long-standing suggestion of a link between fractionaliza- 
tion and high-Tc superconductivity 0, @| be established; 
recently, it has been shown [§] that the topological degen- 
eracy in these systems might be used to protect quantum 
bits and be applicable to the emerging field of quantum 
computing. 

Fractional statistics, as a generalization of the idea 
of quantum statistics based on Berry's phase Q, is a 
sensible idea only in one or two dimensions, where one 
can define a winding number. In ID, the behavior is 
known to occur in spin- 1/2 antiferromagnets [Tot 11 1, 



where exactly solvable models exhibiting this behavior 

exist 0,(11 EI EI El- 

Fractionalization of statistics 
is also known to occur in 2D in the presence of a mag- 
netic field that violates the discrete symmetries of parity 
(P) and time-reversal (T) ; this situation is realized in the 
fractional quantum Hall effect 0, El, EI EI EH El 
(FQHE). Very recently, the fractional statistics of the 
quasiparticle excitations in the FQHE has been ob- 



served experimentally [2l|, |22j. In contrast to the one- 
dimensional case, however, there has been no definite ev- 
idence as to whether fractional statistics occurs in the 
absence of an external field breaking these symmetries. 

In this Letter, we present a spin Hamiltonian for which 
the chiral spin liquid [23, 24 1 (CSL) is the exact ground 
state. The CSL, one of the paradigmatic systems to in- 
troduce the concept of fractional statistics in 2D spin 
systems, is constructed to spontaneously violate the sym- 
metries P and T; this violation is generally associated 
with fractional statistics. The excitations of the liquid — 
spinons, which carry spin 1/2 but no charge, and holons, 
which carry charge but no spin — obey fractional statis- 
tics. In addition, the spinons exhibit quantum-number 
fractionalization and carry only half the spin of the ex- 



citations in conventional magnetically-ordered systems, 
which are spin-1. In many respects, the Hamiltonian 
we present is a generalization of the Haldane-Shastry 
model [HEl (HSM) to 2D, and provides an exact spin 
model in which fractional quantization can be studied. 
A spin Hamiltonian for a 2D system where the ground 
state is a chiral spin state, but not a liquid, has been 
discoverd by Wen et al. [25(. These authors addition- 
aly argue for the plausibility of a CSL ground state in a 
Heisenberg-like model including six-site interactions; the 
model presented in this Letter is precisely of that form. 

The proof presented below shows that the model has 
an exact two-fold topologically degenerate ground state 
for any number of lattice sites N. This is in contrast 
to models, such as the Rokhsar-Kivelson dimer model 
(RKM) , where the topological degeneracy is only realized 
in the thermodynamic limit [2^|. The exact topological 
degeneracy supports the view that the model will increase 
the accessibility of studying aspects of fractional statis- 
tics in 2D on an analytical and exact footing. At present, 
we have numerically verified the results in this Letter 
by exact diagonalization of the Hamiltonian on a 4 x 4 
lattice. The numerical work confirms that the Hamilto- 
nian has two zero-energy ground states and that these 
are the two topologically-degenerate CSL ground states 
introduced below. A detailed discussion of the numer- 
ics will be presented elsewhere [27j |. In the following, wc 
brie! fly review the CSL ground state, present the exact 
parent-Hamiltonian for the state, and show analytically 
that our positive semi-definite Hamiltonian annihilates 
the CSL ground states. 

Ground state. The CSL was originally conceived by 
D.H. Lee as a spin liquid constructed by condensing 
the bosonic spin flip operators on a 2D lattice into a 
FQH liquid at Landau level filling factor v = 1/2. The 
ground state wave function for a circular droplet with 
open boundary conditions, on a square lattice with lat- 
tice constant of length one, is given by [23, E I 
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where is always referred to as the CSL state. The z's 
in the above expression are the complex positions of the 
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up-spins on the lattice: z = x + iy, with x and y integer. 
G(z) = (— is a gauge factor, which ensures 

that (jT|) is a spin singlet. Lattice sites not occupied by 
z's correspond to down-spins. 

For our purposes, it is propitious to choose periodic 
boundary conditions (PBCs) with equal periods L\ = 
L-2 = L, L even, and with N = L 2 sites. Following 
Haldane and Rezayi [HI], the wave function for the CSL 
then takes the form 

2 M 2 

(z 1 ---z M \i>) = l[# 1 (^[Z-Z v }) Jl^jlzj - z k ]) 

M 

■ HG( Zj )e^-M 2 \ (2) 
i=i 

where M = N/2 and $i(w) = -&t(-w) = ■d 1 {w\e~ v ) 
is the odd Jacobi theta function [29(. The zeros for the 
center-of-mass coordinate Z = Z j must lie in the prin- 



cipal region < Rc(Zi) < L, < Im(Zi) < L and satisfy 
Zi + Z2 = L + iL; the freedom to choose Z\ reflects the 
topological degeneracy and yields two linearly indepen- 
dent ground states for the CSL. These states are spin 
singlets, are invariant under lattice translations, and are 
strictly periodic with regard to the PBCs. 

Hamiltonian. The Hamiltonian for which the CSL is 
the exact ground state is defined in terms of vector op- 
erators S~2j as 

H = ^(fl i -%) + -(f2 i -n j -) , (3) 

m 

where the sum extends over all nearest neighbor pairs 
on the square lattice. The vector operators contain one- 
through-three-site interactions and, in terms of two sets 
of coefficients Kjk and [Zy, are defined as 



= E K i» 
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where the prime on the sum indicates i ^ k. The coeffi- 



cients K^k = K(z k 
are given by 
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j) in the first term of 



K(x, y) = hm 
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where the sum over all full lattice translations zq = 
(£ + i m) L guarantees periodicity in the first argument 
of K. The function P(x, y) is given by 

-,(6) 



p ( , jy) =E Co( ^ [ " ]) 

0<\z o -v\<R 



where Co(x) = cos a; 
malization factor: 



Co(£[x-(y-z )}) n(y) 

cosh a; and where n(y) is a nor- 
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0<\z -y\<R 
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chosen such that P(0,y) = 1. The sums in ([6]) and (JT]) 
enforce the periodicity of K in its second argument. 
The coefficients in the second term of (U]) are given by 

u a = t u (z\- z j- z ^)> where 
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In this expression, the function W(z) is the periodic ex- 
tension of 1/z to the torus: 



tw(t*)= lim E — 



R->oo * — ' Z — Zn 
0<\z \<R 



(9) 



The Hamiltonian ([3]) is constructed to be positive semi- 
definite. Therefore, if fli — ttj annihilates the CSL states 
([2]) , these states will be zero-energy ground states of (J3|) . 

Proof. In order to prove that the vector operator 
£li — ttj annihilates the CSL ground state, we first demon- 
strate that the related tensor operator uji — luj annihilates 
it. Here, u>i is a reducible tensor, i.e., a composition of 
tensor components of different ranks, that may be decom- 
posed into irreducible spherical first-rank (vector) and 
third-rank tensors; the operator &i is the vector compo- 
nent of u>i. The operator = u>f — oj^~, where u>f are 
related by a 7r-rotation about the x-axis, will be discussed 
in detail below after constructing the portion of the proof 
that does not depend on its precise form; it is later de- 
fined as uj^ = Ti + Vi with the two operators Tj and Vi 
given in (|13lTT4)) below. The Wigner-Eckart theorem, in 
conjunction with the fact that the ground state defined 
in j2]) is a spin singlet, guarantees that if u>i — uij is a 
destruction operator for the state, then each of its irre- 
ducible tensor components are also destruction operators. 
Therefore, given that the operator — Uj destroys the 
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ground state, it follows that the vector operator fij 
does as well. 



In order to show that the operator ui% 



is a destruc- 



tion operator for the ground state, we first demonstrate 
the following property: 
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The fact that the function on the right-hand side of (|10[) 
is independent of the site-index j ensures that the differ- 
ence of any two operators u>i — u>j is a destruction opera- 
tor. We consider only nearest-neighbor pairs of operators 
in constructing the Hamiltonian in ^ , as this is the sim- 
plest and most local operator that is also translationally 
invariant. Other choices, however, are possible. 

The reducible tensor operators ujj can be decomposed 
into two operators as uij = u>f — wj , where u>f and u>~ 
are related to each other through a 7r-rotation about the 
x-axis that maps S z and S y into — S z and —S y . The oper- 
ator uj~j will be further decomposed as ui^ = Tj + Vj with 
the explicit forms for these operators given in (flTflfl4")) be- 
low. In order to prove (fT0"|) , we will first demonstrate that 



(z 1 - ■■ z M \ip) 



= f{Z) 
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where f(Z) is an odd, periodic function of the center- 
of-mass coordinate Z. Using the relation between ur^ 
and the invariance of the CSL ground state under such 



a rotation, one can show, without specific knowledge of 
the function f(Z), that 



(Zi ■ ■ -z M \ 



■z M \ip) 



= /(W) 



e {zi ■ --z M } 





1 otherwise. 



(12) 



In the above expression, W = w i ' 1S t ne center of mass 
of the down-spins on the lattice, such that {wi} is the 
complement of {z{\. It is straightforward to show, re- 
gardless of the chosen origin, that the sum of the two 
center of mass terms is a full lattice translation: Z+W = 
(£ + im) L. This means that /(W) = —f(Z) and, given 
the definition of uij above, (fT0|) follows from (fTT|) . 

Having developed the remainder of the argument, it 
remains only to demonstrate (fTTj) for the operator ujf to 
prove that ([3]) is the exact parent Hamiltonian for the 
CSL; this last step is the heart of the proof. The oper- 
ator u>j~ is defined in terms of off-diagonal and diagonal 
contributions as ujf = T; 



Vj where 



^ = \ £ K ijk S+ S- (\ + Sf) 



(13) 



(14) 



with the coefficients defined in terms of the functions in 
([SHH]) above. Considering first the off-diagonal term, its 
action on the CSL ground state may be expressed as 



(zi-- - z M | Tj |V0 = \Y. K ^ ZM I S t S k ( \ + S t) W ■ (I 5 ) 
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This is clearly zero if Zj £ {z\ ■ ■ ■ zjsf} giving half of the equality in (fTTj) . Otherwise, acting onto the bra with the 
operator removes the site Zj and replaces it with the site z k . In addition, the matrix clement vanishes if zi ^ {z\ ■ ■ ■ zj^}. 
Dividing by the wave function yields 

(zi ■ ■ ■ z M | Tj \ip) 1 ^ (zi ■ ■ ■ z l ■ ■ ■ z k ■ ■ ■ z M \ip) 



1 v {z\---z l ---z k ---z M \i>) 



{zx ■ ■ ■ z M \i>) ^ Tl (zx- ■ ■ Zi- ■ ■ Zj ■ ■ ■ z M \i>) 
Using the definition in ([5]) above, this may be written as 

(zi ■ --ZM\Tj \lp) _ 1 ^ P(z - Zp, Z t - Zj) (zj ■ ■ ■ Z t ■ ■ ■ Zj + Z ■ ■ ■ Z M \lp) , ,. 

N lr/ i^j z^O 0<\z o -z\<R U ,Y/ 

where the sum over k has been replaced by a sum over z = z k — Zj . As the wave function is periodic in all of its M 
coordinates, the sums on zq and z may be replaced with a sum on x = z — zq that runs over the entire complex plane. 
Additionally, the ratio of wave function coefficients appearing in (|17p has the form 

{zi ... Zi ... Zj + x ...z M \i>) _ G{x)F{x)e -W (18) 



{Zl - ■ ■ Zi - ■ ■ Zj ■ ■ ■ z M \ip) 

t\ MlI Z - Z i\) t£i ^ U [«J - Z *\) 
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where F(x) is an analytic function of x. Being careful to pick up the points excluded by the sum in (1171) and using 
the above definition of F(x), the action of Tj on the CSL ground state may be written as 



(z 1 ---z M \T j 
(zi ■ ■ -z M \ip) 
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(20) 



This is the major step in the proof since the first term may be evaluated with the singlet sum-rule [3Cj |. This sum 
meets the requirements for convergence (3lj that were not satisfied in the original work, due to the exponential fall-off 
of P{x, y) with increasing x. This gives 



E 



P(X,Zi - Zj) 



F(x) G{x) e 



_ f | X |2 _« 
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dx 



P^X, Zi Zj) 



F(x) 



(21) 



Combining this with the second term in Equation l20l gives 
(Zf-z M \Tj \ip) 



(zi ■ ■ ■ z M \tp) 



where the function of the COM coordinate f(Z), which 
first appears in (|10p , may now be written down explicitly: 



(23) 



L 



The [/-function appearing in (|22[) is the one introduced 
in JS]) when defining the Hamiltonian. The equality here 
is a result of the fact that the VF-function, introduced 
in (|9|) , is related to the logarithmic derivative of the odd 
Jacobi theta functions: 



d , . z — z* 
W{z) = — In i (z) + 

UZ 7T 



(24) 



The operator Vj introduced in (|14p . which only gener- 
ates diagonal terms, is chosen to exactly cancel the sec- 
ond term in (|2"2")l . This proves the identity in (TTTj) and 
from here, the arguments at the beginning of the section 
may be traced backwards to show that fi, — Slj annihi- 
lates the CSL state (|2|), and hence that the CSL is an 
exact ground state of ([3]). The fact that the topologi- 
cal degeneracy is exact for any number of lattice sites 
N in this model is due to the fact that Wi — u)j destroys 
the state regardless of the choice of the location of the 
center-of-mass zeroes Zj. 

Conclusion. We have constructed a Hamiltonian that 
singles out the chiral spin liquid state as the exact and, 
apart from the topological two-fold degeneracy for PBCs, 
unique zero-energy ground state. The proof has been nu- 
merically verified on a 4 x 4 lattice. In analogy to the 
HSM in one dimension, this model provides a framework 
to study spinon excitations and their interactions in a 
two-dimensional spin liquid. For example, one may in- 
vestigate whether the spinons in this model are similarly 
free in the sense that they only interact through their 



fractional statistics, and, if so, whether the many spinon 
states can be classified in similar terms [l6j ]. In any event, 
we have promoted the CSL from an intriguing trial wave 
function to the exact ground state of a spin-Hamiltonian, 
and hence accomplished something analogous to the pro- 
motion of Gutzwiller's wave function [3211 to an exact so- 
lution by Haldane and Shastry TH, 13l|. 
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